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Abstract 

Classical gradient systems have a linear relation between rates and driving forces. 
In generalized gradient systems we allow for arbitrary relations derived from general 
non-quadratic dissipation potentials. This paper describes two natural origins for 
these structures. 

A first microscopic origin of generalized gradient structures is given by the theory 
of large-deviation principles. While Markovian diffusion processes lead to classical 
gradient structures, Poissonian jump processes give rise to cosh-type dissipation 
potentials. 

A second origin arises via a new form of convergence, that we call EDP-conver- 
gence. Even when starting with classical gradient systems, where the dissipation 
potential is a quadratic functional of the rate, we may obtain a generalized gradient 
system in the evolutionary T-limit. As examples we treat (i) the limit of a diffusion 
equation having a thin layer of low diffusivity, which leads to a membrane model, 
and (ii) the limit of diffusion over a high barrier, which gives a reaction-diffusion 
system. 


1 Introduction 

We consider evolution equations ii — V(t, u ) that are generated by gradient systems (GS). 
By a GS we understand a triple (X, £,&), where the state space X is a weakly closed 
convex subset of a Banach space containing the states u(t). The functional £ : X —» 
MU{oo} is called energy, but in applications it may be a free energy, a relative entropy, or 
the negative of the entropy. Finally “Jl is the dissipation potential depending on the state 
u and the rate ii such that D ^(u, u ) G X* denotes the dissipation force. The induced 
evolution equation is the force balance 

0 = D A 0l(u(t),u(t)) + D u £(t, u(t)), (1.1) 
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where the symbol D denotes the (partial) Gateaux derivative or the convex sub differential. 
Quite often, we will use the dual dissipation potential ft* that is defined by the Legendre- 
Fenchel transform of ft(w, •). Then, the evolution equation can be rewritten as 

ii(t) =D{0l*(u(t),—D£(u(t))'), (1.2) 

see Section ETT1 for the details. Since ft and ft* are in one-to-one correspondence, we will 
sometimes denote (X, £, 3?) also by (X, £, ft*), in particular if ft* is given explicitly. 

A third equivalent formulation of the gradient flow is given via the energy-dissipation 
principle (EDP), also called De Giorgi’s (ft, ft*) principle, cf. [DMT801 IAGS05] . This 
states that, under suitable technical assumptions, a curve u : [0, T] —* X is a solution of 
(11.11) or (II.2p if and only if it satisfies the energy-dissipation estimate 

£(u(T)) + @(u) < £(u(0)) with @(u) f ft(u(t), u(t)) + ft*(u(t), — D£(u(t))) dt. (1.3) 

Jo 

We call Q) the De Giorgi dissipation functional. 

A GS is called classical if the dissipation potential ft(w, •) is quadratic, i.e. ft(u, u) = 
^(G(u)u,u) for a linear, symmetric, and positive definite operator G(u) : X —» X*. If 
we want to emphasize that a GS is not classical, we call it a generalized GS. The aim of 
this work is to show that generalized GS arise in two natural ways. First, it is shown 
in [MPR.141 IMP* 15] that they appear via large-deviation principles from a microscopic 
X-particlc system for N —» oo, see Section 12.41 for a brief summary of the main result. 
Second, generalized GS occur as suitable multiscale limits of classical GS. 

Obviously, every GS generates exactly one gradient-flow evolution equation by (11.11) or 
(II.2p . but a given evolution equation ii — V(t, u ) may be generated by many GS. If there 
exists at least one such GS, we say that the evolution equation has a gradient structure, 
if we do not want to specify the particular GS. As an elementary example we treat the 
scalar ODE 

p = 1 — 2p with p(t) G [0,1] = Prob({l, 2}), 

which we interpret as the Kolmogorov forward equation for a Markov process (X t )t>o with 
X t e {1, 2}. Obviously, this ODE is generated by the GS ([0,1], £ 2 , ft 2 ) with 

£2 (p) = a (p- 1/2) 2 and %(p,p) = ^p 2 

for any a > 0. Of course, GS that simply differ by a scaling constant such as a > 0 are 
not considered as different. Motivated by a Markovian large-deviation principle, a truly 
different GS is obtained for a > 0 by 

£m v (p) = a(plogp+ (1—p)log(l—p)) and ft^feO = a\/p(l-p) ^*(£/a), 
where the function ^ and its Legendre dual c <o* are given by 

^(v) = 2v arsinh(u/2) — 2\/4-|-u 2 + 4 and ^*(£) = 4(cosh(^/2) — l). (1.4) 

The functions ^ and ^* will play a fundamental role, so we give some elementary relations: 
tf( v ) = y + 0 (v 4 ), ^*(0 = ^ 2 + O(£ 4 ), <r(u) = 2arsinh(u/2), 

y/pqff* (log p— logg) = 2 (y/p - y/q) 2 , y/pq (^f)'(logp- log q) = p - q. 
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Indeed, using the last relation and D£ Mv (p) = oflogp — log(l— p)) we easily find 

P = D € ft M v(p, -D£ Mv (p)) =1-2p. 

Moreover, using = 2sinh(£/2) we see that the evolution takes the exponential 

form 

P = -2y/p(l-p) sinli (-D£ Mv (p))- 

This form is derived and extensively studied in [ BoP 11 . it occurs in mechanics [RRGOO , 
Eqn. (5)] and in chemistry, see the discussion at the end of Section 12.3.21 

The usage of generalized GS is common in the modeling of materials, e.g. for plasticity, 
ferromagnetism, etc., where the nonsmoothness and nonlinearity of the constitutive law 
u n- u ) for the dissipative forces is essential, see Section 12.3.11 and the survey 

|Miel5bj . The mathematical usage of generalized GS in smooth models such as reaction- 
diffusion equations and systems is rather new. One of the remarkable origins of gradient 
structures arises from the interpretation of a macroscopic system as a Kolmogorov forward 
equation 

p = Q* P, where pit) 6 Prob(S'), (1.5) 

for a Markov process (A"(f)) t > 0 on the set S with generator Q. Considering N independent 
particles Xj(t), j — 1,..., N one can define the empirical process p N (t ) = A G 

Prob(S'). For N —* oo the process p N converges to a solution p of (ll.5|h Moreover, 
according to the program in [ AD*111 AD*13) p N satisfies a large-deviation principle that 
gives rise to a rate functional J^(p(-)) = f Q £(p(t), p(t))dt, where £ can be characterized 
explicitly by Q. The main observation in IMPK 1 ll is that £ defines a GS (Prob(S'), £, A) 
via the explicit representation 

£(p,p) = A(p,p) + A*(p, -D£(p)) + D£(p)[p], (1.6) 

whenever Q has a unique steady state n G Prob(S') and Q satisfies the detailed-balance 
condition with respect to 7r (i.e. the Markov process is reversible). We refer to Section 
I2.4l for details, where we also highlight that the arising gradient systems are classical only 
in the case of diffusion processes. In case of jumps, one obtains generalized GS involving 
the function In particular, for p = 1—2p one finds ([0,1], £m v , 3Riv) with a = |. 

We consider the above stochastic approach as a first microscopic origin of GS. The 
second origin involves the concept of evolutionary T-convergence for GS, see the surveys 
|Serlll IMiel5aj for the general ideas. Here we concentrate on convergence results based on 
the EDP, cf. (I1.3|h which is an ideal tool for doing a limit passage for solutions u E : [0, T] —$■ 
X for a family (X, £ e , A e ) of GS depending on a small parameter e. The aim is then to 
derive a limiting GS (X, £ 0 , A 0 ) such that a limit u of the solutions u e is indeed a solution 
for the limiting GS. Our Definition 13.21 introduces the concept of EDP-convergence: A 
family of GS (X, £ e , fk £ ) converges to the GS (X, £ 0 , fk 0 ) in the EDP sense, if the following 


3 

























holds: 


u £ : [0, T] —» X is a 
solution of (X, £ e , lR e ), 
u e (0) —*■ w°, and 
£ E (0,u E (0)) ->• £ 0 (0,u°)<oo 


{ 3u sol. of (X, £ 0 ,1R 0 ) with m(0)=m° 
and a subsequence £k —* 0 : 

Vi 6 ]0,T]: u £k (t) —*■ n(i) and 
£ ^(w fc (*)) ->• £o(«(i)); 


£, 


£0 in X; 


(1.7a) 


(1.7b) 


u E (-) w(-) in L°°([0, T]; X) and 

sup £ e]o,i],te[o,T] &e{u e {t)) < C < 00 


liminf ^ E (« E ) > @o(u). 
6 —^0 


(1.7c) 


When asking only for condition (11.7 a p we speak of pE-convergence, see Definition 13.11 
Note that (ll.7cl) enforces a liminf estimate of De Giorgi’s dissipation functionals Q> £ along 
general functions u £ , not only along the solutions of the GS (X, £ £ , 0l e ). Having this liminf 
estimate, it is easy to pass to the limit in the e-dependent energy-dissipation estimate 
(E3, since the initial energy on the right-hand side is assumed to converge according to 
(ll.7ajh Then, applying the EDP for the limiting GS (X, £ 0 , 0l o ) we see that u is a solution. 

In fact, many approaches to evolutionary T-convergence establish EDP-convergence, 
but do not explicitly state condition (II,7c[) as a main result. E.g. the Sandier-Serfaty 
approach |SaS04 » ISerllj , where the terms / Q T 0i £ dt and / Q T 01* d t are treated separately, 
provides EDP-convergence. Our approach is more general than the latter, since we only 
ask that the sum f Q 0i e d£ + f 0 01* dt behaves well, but not necessarily the individual 
terms. This has two effects: (i) we can allow for general functions u £ , and (ii) it can 
lead to exchanges between the two terms in the limit £ —> 0. Point (i) is important 
to explore outside of the set of solutions and thus providing the full information 
about the GS (X, £ 0 ,fk 0 ), while the set of solutions of the limit equation ii = Vo(t,u ) : = 
D^!ko(u, — D u £o(t, u)) only contains information on Vo- Point (ii) is relevant for another 
important message of this paper. The EDP-limit of classical GS can be a generalized GS. 
This phenomenon is considered as another microscopic origin of generalized GS. 

Here we provide three different examples for point (ii), the first of which is an ODE 
example in Section 13.3.31 while Sections 0] and 0 contain more elaborate examples treating 
the membrane limit of a thin-layer and the limit of diffusion to reaction, respectively. 

For the membrane limit we consider a diffusion equation with a thin layer with very 
small diffusivity. In [Liel2L ILiel3j pE-convergence to the membrane limit was established; 
however EDP-convergence was not studied. We start with the diffusion equation which 
is the gradient flow for the classical GS (Prob(D), £,0l*) with 8(u) = f Q ulog(2u)dx and 
x:(u,o = §/ n a« (x)(d x {;) 2 udx. Using suitable scalings for the diffusion coefficient a £ (x) 
Theorem 14. 1 1 provides EDP-convergence to the generalized GS (Prob(D), £, !Rq) with 

/ 0 p 1 

^(d^fudx + a*V m (° - M 0+ ) ^*(£(0 + )-£(0~)) + J |( d x £) 2 udx , 


where w(0 - ) and w(0 + ) denote the limit of u(x) at x = 0 from the left and from the right, 
respectively. Thus, 01 g involves tf* and is therefore non-quadratic. 

Section [5] follows [ PSVTQl IPSV121 lAM*12] by considering the limit from pure diffusion 
in physical space and along a reaction-path variable y G T cR to a limit of a reaction- 
diffusion system on D. The Fokker-Planck equation reads 

u = mnA x u + r £ d y (d y u + ^ u d y V (y)j , 
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Figure 1.1: For reversible, time-continuous Markov processes the large-deviation principle 
(LPD) of Section l2~il provides a (generalized) gradient structure. This mapping commutes 
with taking the limit e —> 0 and EDP-convergence, respectively. 


where Id is a potential with two global minima y 0 and y\ and one global maximum in- 
between. This equation is generated by the classical GS (Prob(DxT), £ e , Ik*) where £ e is 
the relative entropy and Ik* is the quadratic Wasserstein dissipation potential, see (I5.2|h 
Theorem 15.21 establishes EDP-convergence to a generalized GS (Prob(f2 x {y 0 , y\}), E, R*), 
where R* again involves the non-quadratic function c £*. 

We conclude our introduction by a general and surprising observation. The three main 
models in this work (i.e. the ODE, the membrane, and the reaction-to-diffusion model in 
Sections 13.3.21 [4l and [H respectively) can be seen as Kolmogorov forward equations for 
naturally associated Markov processes. Thus, the large-deviation theory of Section 12.41 
is applicable and provides entropic GS (Prob(S'), £ e , fk e ) for the associated Kolmogorov 
forward equations p £ = Q *p £ for each e e ]0,1[ as well as for £ = 0. The limit for £ = 0 
can be also defined in terms of the classical convergence for Markov processes asking 
p e (t) = e t( ®®p(0) —*■ p(t) = e tJ ®°p(t). Ignoring the (linear) Markovian structure, we can 
also consider EDP-convergence of the induced entropic GS. In all our three examples we 
find the surprising result that the EDP-limit is exactly the entropic GS of the limiting 
Markov process. This means that applying the described large-deviation principle and 
taking the limit £ —> 0 (either on the level of Markov semigroups or as EDP-convergence 
for GS) commute, see Figure [Li] This result appears naturally, if we use representation 
(11.611 of the rate function giving 

S(s) = I + Vtft-DEW) + D£(p)(p]d« = 9{p) + £(p(T)) - £(p(0)). 

Jo 

Hence, the above large-deviation principle exactly encodes the energy-dissipation princi¬ 
ple, and EDP-convergence for the induced entropic GS can be interpreted as T-convergence 
of the rate functionals. 

The question how general this observation about the interchangeability of the suitable 
large-deviation principles and the EDP-convergence is, seems to be challenging, but goes 
beyond the scope of this work. We mention that in [BoP14] similar relations between 
large-deviation principles and evolutionary T-convergence are studied. 

As a final general remark, we emphasize that this paper focuses on the modeling 
aspects of the emergence of generalized GS. Thus, we do not give the full analytical 
details in terms of estimates and convergences in the proper functional spaces, but rather 
highlight the structures and manipulations needed to understand the corresponding limit 
procedures. 
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2 Classical and generalized gradient systems 

We now convert the formal ideas from the introduction into rigorous mathematical state¬ 
ments. We call a triple (X, £,(R) a gradient system (GS), if X is a Banach space, 
£ : [0,T| x X —> Mqq := MU{oo} is a functional (such as the free energy, the negen- 
tropy, etc.), and 01 : X x X —» [0, oo] is a dissipation potential, which means that for all 
g 6 X the functional fR(w, •) : X —)■ M^ is lower semicontinuous, nonnegative, convex, 
and satisfies Ol{u, 0) = 0. In this section, we allow for the case that the energy functional 
depends on the time variable t G [0, T\ to show that the abstract principle is valid in this 
general case. However, for notational convenience we will restrict to the autonomous case 
(i.e. d t E{t,u) = 0) in all other parts. 

We speak of a classical GS, if 0l(u, •) is quadratic, i.e. there exists a symmetric and 
positive definite operator G such that 0l(u,v) = ^(G(u)v,v). However, plasticity requires 
non-quadratic dissipation potentials, e.g. of the form (R(7 t) = cr yie idII ttIIl 1 + l^viscIMI^!, see 
[ Mie03l lMiR15j . In particular, the rate-independent case is based on 0l(u, Xv) = A 0l(u,v) 
for all A > 0, which is incompatible with a quadratic form. If 0l(u, •) is non-quadratic, we 
call (X, £,1R) a generalized GS. 


2.1 Variational principles for gradient systems 

The following proposition from convex analysis shows that there are several completely 
equivalent formulations of the generalized force balance (11.11) . The equivalences of the 
points (ii) to (iv) below are also called Fenchel equivalences, cf. |Fen49j . The essential 
tool is the Legendre-Fenchcl transform T* : X* —y M^ of a convex function : X —> M^ 
defined via 

:=sup{(£,v)-tt(u)|uex}. 

In a reflexive Banach space we have (T*)* = T. 


Proposition 2.1 (Equivalent formulations) Let X ie a reflexive Banach space and 
T : X —> Moo be proper, convex, and lower semicontinuous. Then, for every (6X* and 
every u£X the following five statements are equivalent: 

(i) v G Argmin (T(ta) — (£,w)); (ii) f E iy)\ 

™ex 

(Hi) 'F(u) + tf*(0 = (f,u); 

(iv) v G <9T*(£); (v) £ G Argmin (d / *(^) — 

J7GX* 

Note that the definition of T* immediately implies the Young-Fenchel inequality H /(w) + 
v k*(?7) > ( 7],w) for all w and ip Thus, (iii) expresses an optimality as well. 

Defining the dual dissipation potential 3?* via 01*(u, •) := (fR(w, •))* we can apply these 
equivalences to reformulate (11.11) in the following ways: 

(I) Rayleigh principle |Ray71| 

(RP) u G Argmin (oi(u,v) — (D E(t,u),v)\, 
vgx A ' 


(II) Force balance in X* 


Rayleigh-Biot equation |Ray71 [ fBio55] 


(FB) 0 G du0i(u, it) + D£(t, u) G X*; 
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(III) Power balance in M De Giorgi’s (3?, ft*) formulation [DMT80] 

(PB) ft(u, u) + ft*(w, —D£(f, u )) = — (D£(t, u),u ); 

(IV) Rate equation in X Onsager equation [Ons3l] 

(RE) u e ^ft*(u, — D£(t, u)) 6X; 

(V) Maximum dissipation principle cf. e.g. [HaF08j 

(MDP) D£(t,u) G Arg max ((£, u) — ft* (u, £)). 

£ex* ' ' 


In fact, |Ray71 Eqn. (26)] also includes the kinetic energy ft, which we omit in our ap¬ 
proximation, namely ^(D„ft(u, u)) + D„ft(w, w) + D q £(t,u) = 0. 

Before returning to the general situation, we highlight the three different cases (II)- 
(IV) for the classical viscous dissipation, i.e. H(u,v) = ^(G(u)v,v) and ft*(u,£) = 
|(£, K(m)£) with K(tt) = G(w) _1 . Then, we have 


(FB) G(«)m = -D£(m) (RE) u = -K(u)D£(u) =: -V g £(m) 

(PB) -(G(u)u,u) + -(D£(u), K(m)D£(m)) = — (D£(m), ft), 
where (RE) can be seen as a 11 gradient-flow equation ”, as Vg is the gradient operator. 


2.2 The energy-dissipation principle 

The above formulations can already be understood in a variational sense, since the evo¬ 
lution is expressed by extremizing a functional or by variations or derivatives of the two 
functionals £ and ft. However, for mathematical purposes it is desirable to have formula¬ 
tions in terms of a minimization problem for the whole solution trajectories u : [0, T] —> X. 
One such principle can be derived on the basis of the power balance (PB) by integration 
in time and using the chain rule and finally employing the Young-Fenchel inequality 
ty(w) + T *(rj ) > (r],w), cf. [DMT80] or the survey [MielSaj . This leads to the celebrated 
energy-dissipation principle, also called De Giorgi’s (ft, ft*) principle, see [ AGS05 j for 
some historical remarks. 

Theorem 2.2 (De Giorgi’s energy-dissipation principle) Under suitable technical 
conditions on (X, £, ft) a function u : [0, T\ —> X satisfies (I)-(V) from above for almost 
all t G [0, T] if and only if the energy-dissipation balance (EDB) holds: 

{ £(T,u(T))+ f ft(n, u) + ft*(w, — D£(t, u)) dt 
Jo T 

= £(0,w(0)) + J d t £(t,u(t))dt. 

Under additional technical conditions it is sufficient to have only the upper estimate where 
is replaced by In this case, we speak of the energy-dissipation estimate (EDE). 
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2.3 Examples of generalized gradient structures 

Here we give some examples of generalized gradient structures. First, we discuss dissipa¬ 
tive material models like plasticity or shape-memory materials that form a huge class of 
generalized GS. Second, we treat nonlinear reaction-diffusion systems (RDS), which will 
be closer to the main theme of this paper. The third class of examples concerns reversible 
Markov processes, where the Kolmogorov forward equation has a gradient structure with 
the relative entropy as energy functional. This latter class is so important that it is treated 
in the subsequent Subsection 12.41 

2.3.1 Dissipative material models 

The state of a body C M d , composed of so-called dissipative materials (also called 
standard generalized materials), is given in terms of the elastic displacement u : 12 —y 
and an additional internal variable z : H The latter may describe plastic 

deformations, damage, phase-field variables, magnetization, or other internal states of the 
material. The total stored energy £ depends on u, z, and usually also on a process time 
t G [0, T], i.e. £(f, u, z). As introduced in the theory of standard generalized materials 
in II la.\75 the dissipative forces are given in terms of a (primal) dissipation potential Ik 
that also may include viscoelastic terms: 

fk(u, z ) -k c jiss('2') T -k v j sc (e(u)), 

where e(u) = ^(Vu+(Vu) T ). As before, the corresponding force balance equations (FB) 
are 

0 = Du!k(u, z) + D u £(f, u ,z), 0 G <9ilk(u, z) + D z E(t, u, z) 

with <9ifk(u, i) denoting the set-valued convex subdifferential. 

While the viscoelastic potential !R visc is assumed to be quadratic in many applications, 
the potential Ikdiss for the internal variables z is often supposed to be non-quadratic. E.g. 
in viscoplasticity with yields stress <7 yie id one takes the form 

-^diss(^) J" ^Oyield | % | 1 T v\z\ ^ dx, 

where 6 > 0 is usually taken small, e.g. in 5 = 0.012 in [ ZR*06j . The weak growth of 
order 1+5 is sometimes even replaced by a growth 0{\z\ log |i|) as given by our function 
^ (see e.g. [RRG00 , Eqn. (5)] and [ BoP14] ). 

For later reference, we mention the very simple scalar hysteresis model of a so-called 
play operator. It is given by the generalized GS (M, £ p i ay , fk p i ay ) with 

£piay {t,z) = -z 2 — £(t)z and fk p i ay (i) = r\z\ with r > 0. (2.1) 

It serves as a limit for evolutionary T-convergence in Example 13.31 as well as a large- 
deviation limit in [ BoP 14] . 

2.3.2 Nonlinear reaction-diffusion systems 

We consider concentrations c (t) : O —>- [0, oo[ / of chemical species C\,... Cj that can react 
according to R reactions of mass action type given by a stoichiometric relation 


a\C x + ... + a^C: & ^ + ... + /3JC7, 

JUb 














where r — 1,..., R is the index of the reaction, kr r and A0 are the forward and backward 
reaction coefficients, and the stoichiometric coefficients a\ and are nonnegative integers. 
The reaction-diffusion system (RDS) for the concentrations c = (ci,..., Cj) takes the form 

R 

c = DAc - R(c) with R(c) := ^ ( a r -(3 r ) (2.2) 

r= 1 


and D = diag(5j)j=i ) ... i 7, where Si > 0. With the stoichiometric vectors ct r = (c00 and 
/ 3 r = G Nq we define the monomials in the form c“ := ]/[ i=1 C T■ 

It was shown in [Mielll IMiel3b] that (12.2(1 has a (classical) gradient structure under 
the additional assumption of the detailed-balance condition, which means that 

3 w = (wi)i : Wi > 0 and k { r w“ r = for r = 1, ..., R. (2.3) 


Using the Boltzmann function Ab(^) = zlogz — z + 1 we define the relative entropy 


£(c) 



Xb^x) / WijWi&x, 


which gives rise to the vector of thermodynamic driving forces (also called chemical po¬ 
tentials) n = D£(c) with fii = log(cj/wj). Because of the logarithm laws they satisfy 
the relation at r ■ D£(c) = Y2i=i a i^°s(, c i/ W i) — log(c“ r ) — log(w ar ). Thus, using the 
detailed-balance conditions we obtain the relation 


(a r -f3 r ) • (—D£(c)) = log(^c /3r ) - log(^c Qr ). (2.4) 

To construct the dual dissipation potentials we may choose any scalar, strictly convex 
dual dissipation functional 0 : R. —» M. with 0(0) = 0'(O) = 0 and 0"(O) > 0 and let 


R 


: = / hE fc i v ^i 2 + A3 


2—1 


r =1 


dx 


e c p r 


ki c“ r 


^ ^ 0'(log(/c^c^ r )—log(^c“ r )) 


Using D m 1R(c, /0 = -(div(^QV/Xi)).+^f =1 R^(c)0 , ((a r -/3 r )-/Lt) (a r -(3 r ) and([2“4]), 
we easily see that the nonlinear RDS (12.21) satisfying the detailed-balance condition is 
generated by the pair £ and 30 i.e. c = DAc — R(c) = D Al ik(c, —D£(c)). Thus, we have 
found a family of generalized gradient structures of the nonlinear RDS (12.21) . 

The case of quadratic (R(c; •), i.e., 0(?0 = 70/2 was introduced in [Miellj . leading to 
the logarithmic means H r quadr (c) = A(/0c° r , k^c^) with A (a, 6) = (a— b )/(log a— log b), 
see also [Maalll !ErM12l IMaMlbaj . However, it was already criticized in the 1930s that 
the linear relation c = — K(c)/i (i.e. tR*(c, /z) = A(/z,IK(c)^) is quadratic) arising from 
Onsager’s principle is not suitable for chemical reactions if one wants to model systems 
that are not very close to thermal equilibrium. As a solution Marcclin and de Donder 
introduced exponential dependencies between /z and c, see |Fel72[ Def. 3.3] or CIK'OO. 
Eqn. (11)]. In [GrmlO ] Remark iii on p. 77 gives some historical comments and Eqn. (69) 
explicitly features an exponential dissipation potential 5 involving the function (e^ 2 + 
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e 2 — 2). Since the choice ip(£) = %?*(£) is central for our paper, we give 3?* explicitly 
for this case, viz. 


£*(c ,M) = 


/.IS¥' 


VcJ 2 + 


/? 

E 

r=l 


V / M^r((a r -/3 r ) • At) 


dx. 


(2.5) 


We will see that exactly the same structure, up to a trivial scaling factor 1/2, arises via 
the large-deviation principle described next, see also [ MP*15j . 


2.4 Markov processes, large deviations, and GS 

Here we give a rough sketch of the theory in [ MPR14 ] about gradient structures for 
the Kolmogorov forward equation p = Q *p of Markov processes satisfying a detailed- 
balance condition, which are also called reversible Markov processes, for short. The idea 
that large-deviation principles generate gradient structures goes back to |OnM53j (see 
Eqn. (4-21) therein for a quadratic version of the energy-dissipation principle derived by 
large-deviations, called Boltzmann’s principle). The mathematical theory was developed 
only recently, see [AD*lll IAD*13l IMPR14) . 

In Section 12.4.11 we Erst describe a time-dependent large-deviation principle for general 
Markov processes providing a formula for the rate function J^(p(-)) = f Q £(p(t), p(t)) dt 
and then present the result of [MPR14j . which shows that for reversible Markov processes 
the functional L is induced by an EDP for a GS (Prob(S'), £, 3?). In Sections 12. 4. 2l to 12.4.41 
we then discuss a few applications of the abstract result in Theorem 12.31 

2.4.1 Gradient structures obtained via large deviations 

We consider a compact metric space S and denote by Prob(S') the subset of probability 
Radon measures on S equipped with the narrow convergence defined by duality with 
continuous, bounded functions. The Kolmogorov forward equation p = Q *p describes the 
evolution of the probability laws p(t) of a Markov process (A*) t >o, if the law of A" 0 is given 
by p(0) G Prob(S'). The Markov generator is given as Q acting on functions on S, while 
its dual Q* acts on measures such that f s (Qf)dp = f ? fd(Q*p). 

Considering N independent realizations (xf ))t> 0; i = 1,..., N of the underlying 
Markov process, the measure-valued empirical process p N (t) := e P r ob(S') 

can be defined. Using the law of large numbers the limit A^ —» oo gives p N (t ) p(t), 
which solves the Kolmogorov forward equation p = Q*p, see e.g. |Renl3l . Thm. 2.3.1]. Un¬ 
der suitable assumptions, see |MPR14j , it is shown in [FeK06] that the empirical process 
p N satisfies a large-deviation principle with a rate function J^(p(-)), i.e. 

F(p N (-) « p(-)) ^ e - Ny W\ 

see the above references for the proper definition of The main result is that J? 

has the form y(p(-)) = j n £(p(t), p(t)) d t, where L is the Legendre transform of the 
functional !K(p, •) (i.e. L(p,v) = sup f a £dv — 34(p,^) ) given via the explicit formula 

K(P, «):= Qe«)(«)dp(s). 

Js 

We emphasize the simplicity of this formula and the (separate) linearity in p and in Q. 
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A main observation in [MPR14J is that the deterministic case, which is given by the 
relation 

S(p(-))= [ £(p(t),p(t))dt = 0, 

Jo 

can be interpreted as an energy-dissipation principle if and only if the Markov process is 
reversible, which is the same as asking for the detailed balance condition (cf. (12.311 ) for the 
linear Kolmogorov forward equation p = Q *p. Hence, we now further assume that there 
exists a stationary measure n G Prob(S') which has, without loss of generality, the full set 
S as its support. We say that Q satisfies the detailed balance condition with respect to 7r, 
if 


/ /(Qfi , )dvr= / p(Q/)d?r (2.6) 

Js Js 

for all / and g in the domain if Q. Choosing g = 1, we hnd that Q*7 t = 0, i.e. the detailed 
balance condition implies the stationarity. 

This version of the detailed-balance condition for Markov processes coincides with 
the detailed-balance condition for chemical reactions in (12.311 . Indeed, if the ODE case 
c = —R(c) is linear, i.e. c = Ac with Aw = 0, then (12.311 means A ijWj = ASetting 
S — {1,...,/} and Q = A* gives (12.61) . 

In the sequel we will use the Radon-Nikodym derivative of p with respect to i r denoted 
by / = ^ G L> 0 (S I , 7r) and defined via f B ldp = f B fdn for all Borel sets B C S. 


Theorem 2.3 ([ MPR14 1 Sec. 3]) If the Markov process ( X t )t>o on S is reversible, i.e. 
a stationary measure tt G Prob(S') satisfying the detailed-balance condition (12.61) exists 
for the Kolmogorov forward equation p = Q *p, then the large-deviation rate functional 
f 0 £(p,p)dt has the form of an energy-dissipation principle, namely 


£(p, v) = 3l(p,v) + 3?*(p, — D£(p)) + f D£(p)du, 

Js 

where the gradient structure (Prob(S'), £, !R*) is given by £(p) = | dn and the 

dual dissipation potential 

R*(p,Z)= I ((\/7 e_e )(Q(\/7 e ^))--s/7(Q\/7))d7r, where / = ^- ( 2 - 7 ) 

The cited reference contains not only a full proof, but also specifies under what assump¬ 
tions this implication is in fact an equivalence, i.e. the existence of a gradient structure 
implies the existence of a steady state satisfying the detailed-balance condition. 

Since the arguments and proofs in |MPR14j are quite involved, we highlight here the 
main structures and formal calculations to see that (Prob(S'), £, 3?*) is a GS and that it 
generates the Kolmogorov equation p = Q*p. 

We first observe that Ik* is defined in terms of TC via 

= ^(p,£+^log/) - ‘M i log/^, where ^log/ = D£(p) = ^ !og^- 


Obviously, Ik*(p, •) is convex if and only if J~C (p, •) is convex. The latter is independent 
of the detailed-balance condition and can be established as follows. Consider the Markov 
semigroup P t = e t( ® for t > 0. For fixed p G Prob(S') and t > 0 define 


: = [ e 5 P 4 (e ? )dp= f e ax)H(y) p t (x, dy)p(dx), 
Js JsxS 
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where p t with pt(x, •) G Prob(S') denotes the time-dependent Markov kernel. From the 
convexity of £ i —> e _ £ 0)+£(0 and the nonnegativity of p t and p, we conclude that £ i —> .At(£) 
is convex. Using \{PtP — rj) —> Qrj and ,Ao(£) = 1, we see that 

-A)( 0 ) = Jim 7(^(0 - !) 

t —5-0 t t — ^0 t 

is also convex in £. 

By dehnition we have 01*(p, 0) = 0, and the detailed-balance condition implies the 
time reversibility lR*(p, —£) = 01* (p, £). Since this implies DlR*(p, 0) = 0, convexity gives 
the positivity 01* (p, £) > 0. Thus, OZ* is indeed a dual dissipation potential. 

To derive the induced gradient-flow evolution for the GS (Prob(S'), £, 0Z*) we observe 

D £ S*(p,-D£(p))[, ; ] = D £ M(p,0)h] = / (e°(-r ) )Q(e°) +e°Q(e°,))dp= /,d(Q», 

where we used Q1 = 0. This provides p = DgOZ*(p, —D£(p)) = Q*p, which is the expected 
Kolmogorov forward equation. 


2.4.2 A finite-state Markov process 

We consider the finite state space S = {1,..., 1} such that 

Prob(S') = { p = c = (ci,.. ■, c 7 ) G [0, l} 1 | Y!i=i °i = 1 }• 
The Kolmogorov forward equation is the ODE 

c = Ac with A e M /x/ . 


Note that the Markov generator is given by Qfi nite = A T , and the conditions for a Markov 
generator are 

1 

A ij > 0 for all i 7 ^ j and Vi = 1,..., / : 0 = A ji. 

3 = 1 


We further assume that there is a unique positive steady state tt = w 6 Prob(S') such 
that the detailed-balance condition holds, namely A ijWj = Ajjwy. 

Thus, the induced energy functional is £fmite(c) = | 2i=i c i^B(cj/wj). To calculate 
the dissipation potential we use that Q can be split 


1-1 1 

inite EE Q^ with 

i =1 j=i +1 


:= JTiij I — e?®e l + — e*<g)e J ), 


Wj 


Wi 


where rn^ := AijWj = A jpUi and e k denotes the kill unit vector in M 7 . Using the linearity 
in Qfinite of the formula ( 12 .7[) for !kfi nite we can first calculate 


= J] 

k,l =1 


-' 1/2 e- { ‘QU(^) 1/2 e E ' 


= m. 


ij 


(-) 

. Wk 

CiCi 


WiWj 




(-) 

w k ' 


Wk 


1/2 


(e ?1 + e‘ 


A £» _ 


2 )- 
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Summing these terms and using the function ‘if* we find 


1 7-1 1 

Knite( C ^) = g mi i 

i=l j=i +1 


C{Cj 


WiW 


1/2 


^( 2 ( 6 - 6 )) 


and conclude by Theorem 12.31 that the equation c _ Qfinite C is induced by the GS 
(Prob(S'), £finite; -^finite)- 


2.4.3 Linear reaction-diffusion systems 

We now return to RDS as discussed in Section B.3.2l but now consider only linear reactions 
where all stoichiometric vectors cn r and (3 r are given by unit vectors e* and e J , respectively. 
This means that the reaction is a simple exchange reaction Ci ^ Cj. The linear RDS on 
a bounded smooth domain D C takes the form 

c = OAc + Ac, where D = diag(5j)j = i i ...j with S 3 > 0 (2.8) 

complemented by no-flux boundary conditions. The matrix A is as before. Now, Cj(t , •) £ 
L 1 (f2) is the nonnegative concentration of the chemical species Ci. 

This system can be understood as the Kolmogorov forward equation on the state space 
S = Q x {1,...,/}, where the random variable Y t = (X t ,i(t)) undergoes a Brownian 
motion in D with diffusion constant Sj as long as i(t) = j. At discrete times the particle 
can change its type within { 1 ,...,/}, according to the jump process induced by the 
generator Qfinite = A T , and then continue a Brownian motion with the new diffusion 
constant. The full generator is 


i 

(Q f){x,i) = SiAf(x,i) + ^ A kif{x, k), Vf(x,i) ■ v = 0 on dQ. 

k =1 

We now assume that the linear reaction system satisfies the detailed-balance condition, 
i.e. we assume that there is an equilibrium state w with Wi > 0 and A ijWj = AjiWi for all 
i and j. Then, the steady state 7 r £ Prob(S') is given by the product of the d- dimensional 
Lebesgue measure on and w, up to a suitable normalization factor: 

x 1 
7 r = — d x where Z = Wi vol(f2). 

i —1 

By normalizing w suitably, we may assume Z — 1 subsequently. 

Using the Neumann boundary conditions, it is easy to check that the generator Q 
satisfies the detailed-balance condition (12.6(1 with respect to 7 r. 

Hence, we can apply Theorem l2.3l which provides the large-deviation GS for (12.81) . Note 
that p £ Prob(S') is absolutely continuous with respect to 7 r if and only if p = (q dx) i= i^,,j 
with c = (cj)i= i,...,/ £ L> 0 (f2) / . Moreover, ^ = (ci/wi). =1 f shows that the probability 

density jf- equals the vector of relative concentrations Ci/wi. 

The driving functional £ is the relative entropy up to a factor 1/2, viz. 

£(d) = \ J s XB (^) dn = \ ^(52 MciW/w^Wijdx. (2.9) 
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For calculating the dissipation potential 01 we can take advantage of the linearity in 
Q. In fact, © can be split into / diffusion processes and the reaction part, namely 
Q = EL QL ~t~ Qfinite with. Qfinite above. The corresponding functionals J~Cj for 
diffusion processes take the form 

Hj(p,€) = fy [ (|V£(x,j )| 2 + A£(x,j))dp(x,j). 

The dual dissipation potential 01* is obtained by replacing £(x, j) by £ (x, j)+\ log (-y-^) 
and d p(x,j) = Cj(x)dx, where we also use 

/ Vc j (x)-V£(x,j) + c j (x)A£(x,j)<ix= / c,(x)V£(x, j) ■ z/da = 0. 

Jn Jan 

Subtracting the term at £ = 0, writing £ = (£j)j with £j(x) = £(x, j), and using the result 
for Qfi n i t e from above, we arrive at the formula 

&*(c,£) = I X^lWj(z)| 2 c.,(x)dx 
Jn j=1 

!-l 1 

+ / EE - a/ A ki Ci(x)A ik c k (x) (2(£fc(x)-£i(x))) dx. 

i=l k =2 ^ 

This is the same dual dissipation potential as given in (12.51) . except for the factors | and 2 
outside and inside of ( &*. However, these scaling factors arise since in the large-deviation 
result in Theorem 12.31 a factor | appears in the definition of £(p) (see (12.9ft ). which is 
one-half of the usual relative entropy. 

2.4.4 Large deviations for a membrane model 

We consider a diffusion equation in the interval = ] — 1,1 [, where at x = 0 there is 
a membrane giving rise to a transmission condition. The Kolmogorov forward equation 
takes the form (where ' = d t and ' = d x ) 

P = a ±p" for ± x G ]0,1[, 0 = a±p / (± 1), 

a+p'i 0+ ) = b (p ( 0+ ) - P(0~)) = a~p'(0~). 

The last relation means first that the mass flowing out of ]—1,0 [ has to equal the flow 
into ]0,1 [, and second that this flow is proportional to the difference of the densities. 

The invariant measure is n = |dx, and the Markov generator Q mem b takes the form 

(Qmemb/)(x) = d±f"(x) for ±X6 ]0, l[, 0 = d±f(± 1), 

<*+/'(0 + ) = 6(/(0+)-/(0-)) = o_/'(0-). 

The functional (K mem b takes the form 

^memb(p,0 = [ <*- (£" + (CT) pdx + f d + ((" + (CT) pdx 

J]- 1 , 0 [ 4 ] 0 , 1 [ 

with C'(±l) = o and a + C'(0 + )e c(0+) = 6(e c(0+) - e c(0_) ) = a_C'((T)e c((r) . 


14 







Inserting ( = f + f log(2p) and doing an integration by parts using the nonlinear boundary 
conditions one obtains the dual dissipation potential 


£memb(p,0= / a_(OVdo: + >/p(0-)p(0+) ^*(2(e(0+)-e(0-))) + / a + (f') 2 pdx, 

J]-m J]o,n 

which again features the non-quadratic dissipation function c €*. 


3 Evolutionary r~ convergence 


Following the notions in the survey |Miel5aj we consider families of GS (X, £ e , fR £ ) £e ]o,i[ 
and ask the question whether the solutions u £ for these systems have a limit u for £ —> 0 
and whether u is again a solution to a GS (X, £ 0 , 3?o). Ideally, one might hope that it 
is sufficient for £ e and tk £ to converge in a suitable topology to £o and CR-o, respectively. 
Such results indeed exist and can be found in the surveys (Serlll IMiel5aj . However, 
the aim of this work is to highlight the fact that starting with classical (i.e. quadratic) 
dissipation potentials 3? e we may end up with a limiting dissipation fk 0 that is non¬ 
quadratic. Thus, limits of classical GS may be generalized GS. First such examples were 
given in |Miel21 !MiT12j in the context of plasticity. 


3.1 pE-convergence of gradient systems 

We first recall the general definition of pE-convergence, which is a short name for evolu¬ 
tionary T-convergence with well-prepared initial conditions. Hence, the letter “E” stands 
for both, ‘Evolutionary convergence and ‘E’nergy convergence, while the letter “p” stands 
for welhP’reparedness of the initial conditions, i.e., £ e (0, w e (0)) —> £ 0 (0, u°)<oo. 

Definition 3.1 (pE-convergence of (X, £ e ,0l e )) We say that the generalized gradient 
systems (X, £ £ ,0l £ ) pE-converge to (X, £ 0 ,fk 0 ) ; and write (X,8 £ ,Jl £ ) ^ (X, £ 0 , fko) 


u £ : [0, T] —> X 
is sol. of (X, £ £ , 3L), 
« e (0) — 1 u°, and 
£ e (0,u e (0)) ->■ £q(0, n°)<ooJ 


["3 u sol. of (X,£ 0 ,lko) with u(0)=u° 
and a subsequence £k 0 : 

Vt e ]0, T\. u £k (t) —*■ u(t ) and 
£<*(«<*(*)) £o(«(*))• 


if 


(3.1) 


Here u £ — 1 u means the weak convergence in the Banach space X. We emphasize that 
the notion of pE-convergence asks for convergence of both, the solutions and the energies, 
but not of the dissipation potentials. However, using the EDP and the convergence of the 
energies, we easily obtain convergence of the integrated dissipations, namely 

[ ( u e*!)(w fc , — D£ £fc ( , u , £fc ))dt = £ £k (u Ek (0)) — £ £k (u £k (T)) 

J o 

->• £o(u(0)) — £ 0 (u(T)) = I Jl 0 (u,u)+Jll{u,-D£ 0 (u))dt. 

Jo 

A first systematic study of evolutionary r-convergence relying on gradient structures 
was initiated in Sandier-Serfaty SaSO 1 . see also [Serlll i Mielhaj . In this approach one 
derives sufficient conditions for pE-convergence based on a limiting passage in the EDB 


£ e (u e (T)) + / %(u £ ,u £ ) + Jl*(u £ , -D£ e (u e ))dt = £ £ (u £ (0)). 

Jo 


(3.2) 
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We observe that on the right-hand side we have the initial energy, which converges 
to the desired limit because of the well-prepared initial conditions. Thus, to obtain an 
(EDE) for the limiting process it suffices to show three lirninf estimates for the terms on 
the left-hand side, namely 

liminf £ £ (w £ (T)) > £ 0 (w(T)); 

£—>0 

T nT 

Jl £ (u £ ,u £ )dt > / fR 0 (u,h)df; 

Jo 

liminf / “Jl*(u £ , — D£ £ (w £ )) dt > / (Rq(m, — D£ 0 (m)) df 
e_>0 Jo Jo 

Of course, it is sufficient that these convergences hold only along (a subsequence of) the 
solutions u £ of (EDP). In the following subsection, we will generalize this approach by 
keeping the terms Jl £ and fk* together. 


lim inf 

£—>0 


(3.3a) 

(3.3b) 

(3.3c) 


3.2 EDP-convergence for gradient systems 

Here we define a new notion of evolutionary T-convergence for GS that, on the one hand, is 
more restrictive but, one the other hand, gives a more precise information on the limiting 
dissipation potential ^R 0 . We use the fact that we do not need to have the two convergences 
(13.3b|) and (13. 3 cl) separately. Indeed, it is sufficient that only the integral over the sum 
of the two terms in (13.211 converges. This approach relaxes the sufficient conditions for 
pE-convergence substantially, since in the limit e —>• 0 the different parts of the dissipation 
may be distributed differently. 

In the quadratic case we always have equidistribution Ol £ (u £ ,u £ ) = tR*(u e , —D£ £ (w £ )) 
for solutions u £ . So, if (13.3bj) and (I3.3cl1 hold, we will still have equidistribution in the 
limit, but only along the limit solutions, while the limit functionals need not be quadratic. 
In |Miel2j it is shown that the limit of a classical GS can be a rate-independent system, 
where we always have fkg(u, — D£ 0 (w)) = 0, so there is not even equidistribution along 
solutions. 

We also add a strengthening condition to obtain our notion of “EDP convergence” by 
asking for the convergence not only along solutions, but rather along a suitable general 
class of functions u : [0, T] —> X. For this, we associate with the GS (X, £ £ , lk £ ) £e [ 0il ] 
De Giorgi’s dissipation functionals S> £ , which are defined as 

S> £ {u) := / Jl £ (u, ii) + (R*(w, — D£ e (w)) dt. 

Jo 


Definition 3.2 (EDP convergence) We say that the family (X, £ £ , 3G) £> o converges 
in the EDP sense to (X, £ 0 , fk 0 ), and shortly write (X, £ £ , lk £ ) E -^ P (X, £ 0 , (k 0 ) */ 


(X, £ £ , 3L) £5 (X,£ 0 ,fto), 

p 

£ £ — 1 £ 0 in X, and 

u £ (-)J^u(-) in L°°([0, T]; X) and 1 
sup eG]0il])te[0)T] £ e (« e (i)) < C < oo J 


liminf^ £ (u £ ) > S!q(u). 
£—^0 


(3.4a) 

(3.4b) 

(3.4c) 
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We emphasize that in condition (I3.4cp . the functions u £ are arbitrary and need not be 
solutions of the GS (X, £ e , 0l £ ). From this definition we see that the convergence conditions 
(hi.-ill obviously imply EDP convergence. However, we will study cases, where (13.31) does 
not hold, but we still have EDP convergence. 

We emphasize that EDP-convergence is to be expected whenever one uses the EDP 
principle for establishing pE-convergence. Indeed, from general arguments one presumes 
that De Giorgi’s dissipation functional f$ e has (after extraction of a subsequence) a T-limit 
in the form 

S>q(u) = / M 0 (u(t),u(t))dt. 

Jo 

From the lower semicontinuity of T-limits, one expects that M(u, •) is convex. Hence, one 
can define via 

3£m(-u, v ) := M 0 (u, v ) — M 0 (u, 0) 

and hope that it is a dissipation potential. For this, one needs to show (i) the positivity 
3 £m(-u, v) > 0 and (ii) M 0 (u, 0) > —D£ 0 (w)). Often, the positivity (i) follows simply 

from the evenness M 0 (u, —v) = M 0 (m, v) and convexity. Moreover, if it is possible to show 
M 0 (u, u) > — (D£ 0 (w),u) (this holds for e > 0 in the form dl £ (u, u)+1R*(m, — D£ e (w)) > 
— (D£ e (w), u)), then we hnd (ii) via the estimate 

— D£ 0 (w)) = sup ({— D£ 0 (w), v) — M 0 (w, v) + M 0 (w, O)') < M 0 (w, 0). 

Thus, we arrive at the desired EDE £ 0 (u(T)) + + 3^^ d t < £ 0 (u(0)) and EDP- 

convergence to (X, £ 0 , CRjvt) is established. 

It would be interesting to study more generally the relations between pE and EDP- 
convergence. Obviously, showing the lirninf estimate for is the major step in establish¬ 
ing pE-convergence. Hence, it seems redundant to ask for the pE-convergence explicitly, 
yet it is not obvious under what additional condition (e.g. the validity of a suitable chain 
rule) we really can deduce the pE-convergence from the lirninf estimate for @ e . 

We end this section with two examples concerning EDP-convergence. Example 13.31 
shows that the model discussed in jMie!2j satisfies pE-convergence but not EDP-conver¬ 
gence. Example 13.41 emphasizes the fact that pE and EDP-convergence are not properties 
of an evolution equation u = V £ (u) but of a GS (X, £ e , 3L). Indeed, for a given equation 
one may have different gradient structures leading to different limits in the EDP sense, 
which in turn generate different limit evolutions. 

Example 3.3 (pE-convergence without EDP-convergence) We consider the wig- 
gly-energy model introduced in l ACJ96 1. It is given via the time-dependent GS (M, £ e , H £ ) 
with 

£ e (t, u ) = -u 2 — i(t)u + re sin (u/e) and 3£ e (w) = - u 2 . 

For sufficiently smooth loading curves i : [0, T] —> M it was shown in lMie!2l Tlim.3.2] 
that the GS (K, £ e , 3£ e ) pE-converge to the generalized GS (M, £ p i ay , 3^ p i ay ) defined in (12.1 [) . 
Obviously, we have the uniform convergence £ e —* £ p i ay , while ( 3> £ converges to a limit 31 0 
that cannot be written in terms of Up i ay + 3?* lay , see lMie!2\ Prop. 3.1] 

Example 3.4 (Different limit equations) Here, we provide an example of an evolu¬ 
tion equation ii = V £ {u) with two different gradient structures. Both gradient structures 
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have an evolutionary T-limit in the EDP sense, and the surprising fact is that the gener¬ 
ated limit evolutions are different. Thus, EDP-convergence and pE-convergence are not 
properties of the family of evolution equations u = V e (u), hut of the chosen gradient 
structures. 

Consider Q = [0,1] and let X = M> 0 (O), the set of nonnegative finite Radon measures. 
Moreover, consider a continuous periodic function A : M. —)■ ]0, oof such that 

0 < a min := min ye]R A(y) < max yGM A (y) =: a max < oo. 

With a £ {x) = A(x/e) we dehne the simple PDE 


u(t,x ) = —a £ {x)u{t,x), m(0,x) = uq(x ) > 0 for x G 

We introduce two different gradient structures (X, £ e , lk e ) and (X, £ e , lR e ) via 

£ e (u)= [ a £ du, T e (u,£)= [ ^ 2 du, E e (u) = [ — du, T e (u,£) = [ u. 

Jo. Jo ^ Jo a e Jo ^ 

It is shown in !Miel5d Cor. 3.8] that these GS converge in the EDP sense to the limit 
systems (X, £ 0 , lk 0 ) and (X, £ 0 ,lko), respectively, where 


£q(^ 0 / ^min d'U, 




du, R* 0 (u,£) 



U. 

2 


In particular, the limit evolution for the hrst is ii = —a m i n u, while it is u = —a max u for the 
second. This is not a contradiction, but has its origin in the well-preparedness condition 
for the initial data. No sequence can be well-prepared for both systems, i.e. if u £ —*■ u and 
Se( M e) £ 0 {u), then we have £ 0 (w) ^ lim inf e ^. 0 E £ {u £ ), and vice versa. 


3.3 EDP-convergence for an ODE example 

We discuss a very simple example of a discrete Markov process with state space S = 
{1, 2, 3}. The jump rates are such that in the limit £ —> 0 the particles never stay in the 
state 2. Thus, the limiting Markov process has the state space {1, 3} only, see Figure ETT1 
We will start with three different GS, namely (i) the quadratic one, where both £ e and 
H £ are quadratic, (ii) the entropic one with classical Ik*, and (iii) the entropic one with 
the dual dissipation potential defined in terms of If*. The interesting point is that in the 
cases (i) and (iii) the limiting GS obtained via EDP-convergence will still be in the same 
modeling class. However, in case (ii) we will lose the classical GS and obtain a generalized 
GS that cannot be described via ff*. 

We consider the Kolmogorov forward equation (here an ODE) of a Markov process on 
the state space S = {1, 2, 3} given by 

/ -1 l/e 0 \ 

11 = ri+.n I 1 -2/e 1 u, X := Prob({l,2,3}). (3.5) 

V 0 l/e -1 / 

The unique equilibrium w £ = i^{l,£, 1) T satisfies the detailed-balance condition (I2.6|h 
The limit dynamics is easily obtained by setting U 2 = sr giving 

iii = (2+e)(r-u 1 ), ef = (2+£)(u 1 -2r+u 3 ), ii 3 = (2+e)(r-u 3 ). 
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1 



Figure 3.1: Left: Three-state Markov process with high rate of leaving state 2. Right: 
The limit for e —>■ 0 gives a two-state Markov process. 

Thus, in the limit £ —* 0 we find r = (ui+u 3 )/2 and 


hi = w 3 — Mi, 0 = U\ — 2r + w 3 , m 3 = «! — n 3 . 


More precisely, if the initial condition satisfies u £ (0) —* ( p 0 , 0,1— p 0 ), then for all t > 0 we 
have 

u £ {t) —)■ (p(t), 0,1— p(t)) T where p(t) = 1 — 2p(t), p( 0) = p 0 . 

We will study the limit e —> 0 in several gradient structures. For general strictly con¬ 
vex and superlinear functions 0 and ip we consider the GS (X,S £ ,0Z £ ) with £ e (u) = 
YPl-\ w i ( P{ u i/ w !) (recall w e = l,e, 1) T ) and 


2 

^*(+0 = 

i=i 


where a}(w) := 


“i+i/^I+i - V” 1 ! 


Using the fact that v — u satisfies V 2 = —v\ — u 3 we obtain the primal dissipation potential 
0l £ in the form 

%e(u,v) = aUu)'ip(-^j) +a|(u)0(^y). (3.6) 

Indeed, in general 0l £ is an inf-convolution, since 01* is a sum over two terms, ffowever, 
here we can eliminate V 2 and argue as follows: 


%(u, v ) = sup{ £ • v - T e {u, f) | 6 + 6 + 6 = 0 } 

= sup{ (6-6)^i + (6-6)Us - K( u > 016+6 + 6 = 0}, 

which gives the desired result, since 01* only depends on 6 — 6 and 6 — 6- 

We now state the result on EDP-convergence of (X, £ £ ,0i £ ) to (X, £ 0 ,CR 0 ). Since the 
limiting GS can be described on the much smaller set Prob({l,3}), which we identify 
with Y := [0,1], we can formulate the limit GS in terms of the reduced GS ([0,1], E. R). 


Theorem 3.5 We have (X, £ £ , ft £ ) E -= P (X, £ 0 , $o) 


£ 0 (m) 

0Z o (u, u) 


E(p) for u = (p, 0,1—p) T , 
oo otherwise ; 


with E(p) 


\<t>0p) + §0(2—2 p) 


R(p,p) for(u,u) = (p,0,l-p,p,0,-p) T , 
oo otherwise , 
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where R is given via its Legendre dual 


R *(p,v) ■= a (p) + SU P \ in f (a(P,z)ip*(v-r)+a(l-p,z)'il}*(T)) - Z(p,z) \ , 
z> 0 v / J 

where S (p, z) :—a(p,z)'ift*[(ft'(2p)—(ft , (z)) + a(l— p, z) ip* (ft(2—2p)— (ft'(z)), 

a(p,r) : = ( 2 p—r)/{ (ip*)' {eft'(2p)—(ft'{r)') }, and 

cr(p) := inf{ E(p, z) \ z > 0 }. 


In particular, this implies that the limiting ODE p = 1 — 2p is induced by the reduced 
generalized GS ([0,1], E.R), i.e. p = 1 — 2p = D ;/ R(p, —DE(p)). The above theorem 
follows directly from the next proposition and the general theory described in Section 13.21 
For the energy-dissipation principle we consider De Giorgi’s dissipation functional 


: (u) = / M £ (u(t), u(t)) dt with M £ (u, v) = 3l e (u, v) + CR*(it, — D£ e (u)). 


r r r 

Proposition 3.6 We have the T-limits £ e —> £ 0 , M £ —> M 0 , and < 2> e with 


M > x _ / R(p,^)+R*(p,0 / ( 2 P)-0 / ( 2 - 2 P)) for i u , v ) = ip, °, 1-P, G 0, **v) 

ivi 0 [u,v)-< oo otherwise- 


T 


and S> 0 (u) = M 0 (u(t),u(t))dt, where R is given as in Theorem \3.5\ . 


r 

Proof: The convergence £ e — )■ £o follows easily from the explicit form £ ; (u) = 

wf(f)(ui/wi), the convergence w 2 —> 0 , and the superlinearity of (ft. 

To simplify the T-limit of < 2> e we introduce the scaling u 2 = er and use that ii 2 = er 
does not explicitly appear in M e , see tR e in (13.61) . However, the relation u\ + u 2 + u 3 = 1 
now takes the e-dependent form ui+er+u^ = 1. Moreover, defining M e (ui, r, w 3 , hi, M 3 ) = 
M £ (ui,er,u 3 ,u) shows that M e is continuous in e G [0,1] and all the arguments. 

Hence, when considering a sequence of functions u £ : [0, T\ —> X with u e —^ u 0 and 
E e (u £ (t)) <C< 00 we find uo(t) = (p(t), 0,1— p{t)) T and 


lim inf $>Ju e ) > 

e —>0 V ’ ~ 


m(p(t), p(t)) dt with 


I V 

m(p,v) := inf <a(p, z)ift(— --) + a(l— p, z)ift( 


z> 0 


a(p, z) 


'a(l ~P,z) 


)+E(p,z)}. 


Note that we have no control over the variable r(t) = u 2 (t)/e in the limit e —» 0. So, we 
simply minimize over all possible values z G [ 0 , 00 [, which certainly provides a good lower 

p 

bound. Moreover, recovery sequences for the convergence S> £ — 1 can be obtained in 

the form u £ {t) = (p(t), 0 , 1 — p{t)) + e((t)(— p(t), l,p(t) — l), where ((t) is the minimizer in 
the definition of m(p(t) , p{t)). 

Obviously, the definition of a gives the relation a(p) = m(p, 0). Thus, we have derived 
the reduced dissipation potential R in the form R(p, v) = m(p, v)—a(p). Doing a Legendre 
transform with respect to v we obtain the form of R* given in Theorem 13.51 since the sum 
turns into an inf-convolution. ■ 


Next we consider three different choices for (ft and ift. 
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3.3.1 Quadratic energy and dissipation 

First, we consider the case 

4>{ r ) = -r 2 and - 0 ( v ) = -r 2 , 

which gives E(p) = p 2 + (1—p) 2 and a(p,r) = 1 and simplihes all expressions considerably: 

Z(p,z) = \i 2 : p-z) 2 + ^(2—2p—z) 2 , a(p) = (1-2 p) 2 , R(p,r/) = ^p 2 . 

Here, the crucial point in the definition of R* is that the inf-convolution involving 
%b* and r does not involve any dependence on z, so that the term a(p) exactly cancels 
sup , >0 — £(p, z), which is generally not the case. Thus, the limiting GS is ([0,1], E. R) 
where R (p, u) = v 2 is quadratic, and ([0,1]. E. R) is again a classical GS. 

3.3.2 Entropic energy and ^-type dissipation 

Next, we consider the case of the Boltzmann entropy and the dissipation defined in terms 
of i/j = which coincides with Section [TT except for the trivial scaling factor 2: 

(f)(r) — A B (r) = r logr — r + 1 and -0*(£) = ^*(£) = 4(cosh(£/2) — l). 

This gives the reduced energy functional E(p) = |Ab(2p) + 1Ab(2—2 p) and a(p, z ) = 
i/2 pz. In the latter expression and in the dehnition of £ we proht from the interaction 
of Ag(r) = logr and the exponential form of %f*, viz. 

£(p, z) = 2(y/2 p— \fz) 2 + 2{\j2—2p — ^fz ) 2 = 4 — 4 b v \fz + 4z with b p = y/2p+ \j2—2p. 
Minimizing in z > 0 we arrive at 

crip) = 4 - h l = 2 - 4^(1— P) = 2 {y/p ~ V^pf- 
For calculating R* we hrst observe, for a, b > 0, the formula 

inf (a^*(r) + Mf*(£-r) N ) = 4\ (a+b) 2 + ^V*(£) - 4(a+6), 

which follows by writing the left-hand side via ^*(£—r) = 2 e^/x+ 2 e^^x —4, where x = e T , 
and minimizing in x > 0. With a = a(p, z) and b = a(l—p, z) we find 

R *(p, p) = a(p) + sup (4 VzJb 2 + a/p( 1 -p)^*(? 7 ) - \b v yfz - £(p, z)) 

= a(p) + sup bl + v^l-p)^*^) - 4 - 4z) 

z>o 1 v 1 

= <r(p) - 4 + b 2 p + y/p(l-p) = y/p(l-p) ^*(p). 

We emphasize that in this minimization with respect to z it is crucial to keep the terms 
involving the dual dissipation potential ( ^ > *{p) and the term £ together. 

We observe that the resulting gradient structure is again the structure, which is ob¬ 
tained from the large-deviation principle of Section 12.41 This confirms the statement that 
gradient structures obtained from the large-deviation theory are very stable against taking 
further limits in the sense of EDP-convergence, see Figure [1J] 
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3.3.3 Entropic energy and quadratic dissipation 


In [Maalll IErM121 IMiel3al IMaM15al IMaM15b] the relative entropy was used for the 
energy functional £ and a quadratic dissipation leading to a classical gradient system: 

0 (r) = Ab(t) = rlogr — r + 1 and = -£ 2 . 

We obtain the same limit energy E (p) = |Ab(2 p) + |Ab(2—2 p) as in the previous case, 
but the functions a,j(u) are quite different as they involve the logarithmic mean A(r, s) = 
log r- log s * I n deed we have a(p, z ) = A(2 p, z ) = lu „.^~^- . We further obtain the functions 

^ (i^p—z) ( log(2p) - log z) + (2-2p-z)(log(2-2p) - logs)) 


and have no explicit formula for cr(p) = inf 2> o E(p, z). In the dehnition of R* we can do 
the inf-convolution explicitly, since ip* is quadratic, so we find the formula 


R* {Pi v) 


a(p) + sup 

z>0 


a(p, z)a(l—p, z) 
2(a(p, z)+a(l-p, z)) n 



We claim that the growth of R*(p, rf) is no longer quadratic, but exponential. For this 
we insert z = e br ' for rf 1 for some b G ]0,1/2[ into the supremum to obtain a lower 
bound. From T,(p,z) ~ zlogz and a(p,z) ~ z/logz for z —* oo we hnd the asymptotic 
lower bound 

R*(P,V) ~ (-£ - 

Hence, we see that the growth is at least as e b ^ for all b G ]0,1/2[. Moreover, we expect 
that the function R*(p, rf) does not have a product structure 6 (p) v h(?/) any more. 

Thus, we see that the classical gradient structure for the relative entropy is not stable 
under EDP, in general. Nevertheless, in ICli. M 13. lDiL14l Ala Al I da evolutionary T-limits 
between discrete Markov processes and continuous Fokker-Planck equation are studied, 
where the classical gradient structure survives. 


4 The membrane as a thin-layer limit 

In our first major application of the EDP-convergence as a microscopic origin of general¬ 
ized GS, we follow [Liel21lLiel3j and consider a one-dimensional diffusion equation with a 
thin layer of very small diffusivity. Assuming that the diffusion coefficient and the width of 
the layer scale in the proper way, we will arrive at a membrane model in the limit. While 
the limit passage of the linear diffusion problem to the linear transmission problem at the 
membrane can be done directly (or with the quadratic gradient structure, see (Liel3| ). 
we prefer to do the somewhat more elaborate EDP-limit using the GS with the relative 
entropy as energy functional and the classical dissipation potential of Wasserstein type. 
This case was already studied in [Li el 2 i Sec. 3.2] in a more special setting and without 
explicitly calculating fRg. 

We start from the equation 

u = (a £ (x) (V + in Q := ]— 1 , 1 [, d x u(t, ± 1 ) + u(t, ± 1 )V^(± 1 ) = 0 , (4.1) 
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where ' = d t and ' = d x . By our choice of the boundary conditions, the total mass 
J n u(t,x) dx = 1 is conserved, thus we can interpret the equation as the Fokker-Planck 
equation of a Markov process. Defining the equilibrium density 

1 r 1 

w £ {x) = — e ~ v ^ x) with Z £ = J e~ v ^dx, 


we have the GS (Prob(il), £, fk*) with 

£ e (w) = / Ab(u(x)/w £ (x))w £ (x) dx and fk*(w,£) = - / a £ (x)u(x)£'(x) 2 dx. 

Jn 2 J n 

The nontrivial behavior happens in the thin layer given by the small interval [0, e]. In 
particular, we allow a £ and V £ to depend non-trivially on x: We assume that there are func¬ 
tions a*, a+, V*, V + G C 1 ([0,1]) and a_,V4 G C 1 ([— 1,0]) such that a*(x), a+(x), a_(— x) > 
a > 0 for all x G [0,1], 14(0) = 14(0), 14(1) = V+(0), and 


a £ (x) 


a + (x) for x > e, 

ea^x/e) for £ G [0,e], and V e (x) 
a-(x) for x < 0, 


V + (x+e) for x > e, 
vXx/e) for x G [0, e], (4.2) 

V-(x) for x < 0. 


Here 14 is constructed to be continuous on = [—1,1], while a £ has jumps for x G {0,e}. 

The pE-convergence result established in |Liel2l Sec. 3] states that the limiting system 
is given as a membrane problem, where the thin layer is replaced by a transmission 
condition. The interesting point is that the EDP-convergence reveals that the limiting 
GS is no longer classical but involves for the jump of the driving forces at the membrane. 
For passing to the limit we note that the function w £ converges pointwise to the limit 


wo (4) 


J_ e bW for x > 0, 

_i_ e - y -( x ) f or x < o, 


with Z„= I e- v -^dx + J e _K| - w drc. 


(4.3) 


which may be discontinuous at x — 0, but has well-defined limits u>o(0 - ) and ta 0 (0 + ) from 
the left and from the right, respectively. This limit is totally independent of the potential 
14 inside the layer. The influence on the layer potential 14 and the layer diffusion profile 
a* will only survive in one coefficient 4*. 


Theorem 4.1 (Membrane limit) (Prob(f2), £ e , %*) E —^ (Prob(f2), £q, 3?q)> where 


£o (u) 

X*(u,0 


AB(u/w 0 )iv 0 dx and 


a_ 




(^') 2 wda; + 


']o,i[ 


y(< i'fudx 


+ 


yl * \J Z (S : j^o+ T^*(^(° + )-^(° )) where A, 



Z 0 e v *(yl 

a *(y) 



-l 


Before we go into the details of the proof, some comments are in order. First, we 
emphasize that the constant Zq in the definition of the coupling coefficient A * is not 
related to 14, but only depends on V±, see (14.31) . Hence, for a large barrier 14 the 
transmission coefficient H* becomes indeed small. 
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Second, the limiting equation is a PDE in the subdomains = ] — 1, 0[ and G + = 
]0,1[ coupled by a transmission condition. It can be obtained easily by considering test 
functions £ G H 1 (Q_) x H 1 (h2 + ) in the weak form f n iif dx = — D£ 0 (w)) [£]. Using 

the fact that £ may have a jump at x = 0, the transmission conditions arise via the 
boundary terms when integrating by parts. We arrive at 

u = (awo(u/w 0 ))'^ for x G U G + . 

0 = a + (0)n;o(0 + )(«/^o)) / (0 + ) - - ^y), 

0 = a_(0)w 0 (0-)(n/n;o)) , (0") - - ^Fj), 

0 = a±(x)wo(x) (u/w o)) (x) at x G {—1,1}. 

We refer to [ G1M13 ] for a similar derivation of more general nonlinear transmission con¬ 
ditions and active interface conditions using gradient structures. 

Finally, we remark that the primal dissipation potential can be written using the 
integration operator I[u](x) := ff l u(y) dy = — J u(y)dy, where the last relation follows 

from /]) j ii dy = 0, which in turn is due to u(t) G Prob(G). Noting that the functions £ 
may have a jump at x — 0, one has the identity 

[ £hdx = - f /[«]£' dx- 7[h](0)(£(0 + )-£(0-)) - [ /[«]£' dx. 

J -1 J -i Jo 


Here J[m]( 0) is the flux through the membrane, which is thermodynamically conjugate to 
the jump £(0 + )-£(0“) in the driving forces. With this and /[«](—1) = 0 = /[w](l) the 
evaluation of the Legendre transform for IRg yields the primal dissipation potential 


X 0 (u,u) = 


'-i 


m 2 

2 a_n 


dx + 


- 1 m 


2a i u 


dx + A„ 


m(0 )u(0+) 
Ji) 0 (0 _ )lDo(0+) 




J[*](0) 


A„ 


u (0 )u(0 + ) 
wo (0 _ )mjo(0+) 


(4.4) 


Proof of Theorem f.l: We first observe that £ e —>■* £ using [AGS05. Lem. 9.4.2], Moreover, 
pE-convergence was established in {Liel21 Sec. 3.2], It remains to establish the liminf 
estimate fl3.4d) . where De Giorgi’s dissipation functional Qt £ takes the explicit form 


:(u) := 


in 


2 a e u 


/ H 2 + ^(( 1 og©)') dxdt 


Step 1. Blow up = transformation from @ e to Q> £ : To study the T-limit J&q of Q> £ we 
blow up the thin layer such that its transformed thickness becomes of order one. For this 
we use Y e : [—1,1] —> [—1, 2] and its inverse X £ = Yf 1 : 


Y £ {x) 


x for x < 0, ( y for y < 0, 

i±^x for x G [0,e], and X £ (y) = < ^ y for y G [0,1+e], 
x+1 for x > £• [ y— 1 for y > 1+e. 


For u : [0, T] x G —> M. and y G G := ]—1, 2[ we dehne the functions 


U e (t,y) = u(t,X e (y)), W £ (y ) = w £ (X £ (y)), A e (y) 


a e (X e (y)) 

X'(y) 
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and the functionals Q) e via @ e {u) = S> £ (U £ ) and find 


%(U) := / T ^^fjSP] 2 + m~((\og{UIW,))') 2 Ay&t, 

where I e [U](Y) = f\ U{y)X' e (y)dy = - U(y)X' E (y) dr?. 

Following the arguments in (Liel2l Sec. 3.2] it is not difficult to establish the F- 
convergence of 3> e to HZq, where the latter is given in the form 

%(U) := [ + ^((\og(UIW))’) 2 AyAt 

with I 0 [U](y) j U (77) M (77) dr) — — J U(r))M(r))drj, 

_ _ f (a-(y),e~ v -M/Z 0 ,l) for y < 0, 

where (A(y), W(y), M(y)) := < (a*(y), e~ v *^/Z 0 , 0) for y e [0,1], 

( ( a+(y),e- v +( y) /Z 0 ,1) for y > 1. 


Step 2. Minimization over the rescaled layer: The main structure in this limit model 
is that % does not depend on U(t, -) |] 0?1 [, since M(y) = 0 for y e [0,1]. Moreover, 
on this interval J 0 [17(f, •)] is constant, namely y u (t) := J]°j U(jj) dr/ = / 0 [f/(f)](?/) for 
y G [0,1]. Thus, given the value y u (t) we can obtain the optimal profile of 17(f)|[ 0 ,i] from 
the boundary values U(t,Q) and U(t,l) and minimizing the functional S(y u (t)'r) given 
via ^ 

- f 1 rv 2 ATI / _ ,\ 2 

S{a,U) := J ^ + —((log (U/W))) Ay. (4.5) 

Now, Proposition I A. 21 in Appendix I A. II provides the explicit formula 

G(a, u 0 , «i) := min{ S(a, 17) | U > 0, 17(0) = w 0 , U{l)—U\} 



with w- = W( 0) = w 0 (0 ), w + = W{ 1) = w 0 (0 + ), and A* = { l/(A(y)W(y)) dy) 1 
Inserting the definitions of A and W gives exactly the formula for A* in the theorem. 
Thus, we have constructed a simpler functional which is given by 


%(U) := 


0 L 7]-l,0[U]l,2[ '2A17 


1 -mu ?+47 


d y 


+G[Io[U](0),U(t, 0),U(t, 1) 


df, 


satisfies the lower bound S>o(U) > S>q{U) for all 17, and has the important property that 
it does not depend on 17|[ 0 ,t]x]o,i[- 

<Step 3. Relation between and 3> 0 : Using the special form of dig and stated in 
'1 1. II) and the theorem, we dehne the limiting dissipation functional @o(u) = f Q T JZo(u, u) + 
—D£ 0 (w)) df. By construction and the special form of G the functional is closely 
related to SZq in the following way. For any function U : [0, T] x [—1, 2] — y M, we may define 
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u : [0, T] x [—1,1] by u(t,x ) = U(t,Y 0 (x)), where Y 0 (x) = x for x < 0 and Y 0 (x) = x+1 
otherwise. Then, we have S>q(u) = 3>o(U) < 3>o(U). 

Moreover, for any function u one can construct an optimal U as follows. We split u at 
x = 0, the right part is shifted by 1 to the right, and the minimizer of U(t , •) G H 1 ([0,1]) 
of 9(/[h(t)](0), u(t, 0 _ ), u(t, 0 + )) is inserted into the gap. Then, 3> 0 (u) = S>q{U). 

Step 4- The liminf estimate (13.deft : To establish the fundamental liminf estimate 
we consider, w.l.o.g., sequences u £ satisfying 1/R < u £ < R for some large R > 1. 
In particular, by minimum and maximum principles these bounds can be expected for 
solutions of Defining U e (t,y ) = u e (t,X e (y )) we again have U £ (t,y ) G [1/R, R]. 

Thus, we find a subsequence £& —> 0 such that 

u £ u 0 in L 2 ([0, T] x Q) and U £ U 0 in L 2 (0, T; L 2 (D)). 

Moreover, we have uo(t,x ) = Uo(t,Y 0 (x)). Now, using $> £ (u E ) = S> £ (U £ ) we arrive at the 
desired liminf estimate 

liminf @ £ (u £ ) = liminf @ £ (U £ ) > § 0 (U 0 ) > ^o(^o) = %{u 0 ). (4.6) 

£->■0 £—^0 

This concludes the proof of Theorem 14.11 ■ 

We conclude this section by observing that the EDP-limit of the thin-layer diffusion 
system given by the classical GS (Prob(O), £ e ,(R*) is a the generalized GS for the mem¬ 
brane problem. For £ > 0 and for £ = 0 the gradient structures are exactly the ones 
obtained from the large-deviation principle, see Section 12.4.41 Hence, we again found 
an instance where the diagram in Figure 11.11 commutes, that means that applying the 
large-deviation principle can be interchanged with taking the EDP-limit £ —)■ 0. 

5 From diffusion to reaction 

In our second major application of EDP-convergence as a microscopic origin of generalized 
GS, we continue the work in }PSV10| PSV12 . AM* 12] which show that linear reactions can 
be obtained as limits of diffusion for a suitably scaled energy barrier. In jPSVHk PSV12] 
the method relies on a quadratic energy functional and a classical gradient structure. In 
[AM* 12] the pE-convergence for the entropic GS is shown, but only diffusion along the 
reaction path is allowed. In fact, the result therein gives EDP-convergence, if one takes 
the addition in |MPR14l Prop. 4,4] into account. 

Here we generalize the latter work by also allowing diffusion in a physical space D, 
such that the resulting limit equation will be a (linear) reaction-diffusion system. Our 
physical domain O C R d is bounded and has a Lipschitz boundary. For the reaction path 
we choose T = [0,7] C 1 and define the cylinder Q = D x T. (Indeed, T could by any 
bounded or unbounded interval.) 

For densities u G L 1 (Q) the integral f D J'J/ udydx denotes the number of particles per 
unit volume that are in the subdomain D cO and have a reaction state y G [yo,yi] C T. 
The evolution of the density u is driven by diffusion in the ^-direction with diffusion 
constant rriQ > 0 and a much faster diffusion in the y-direction with diffusion constant 
t £ 1 to allow the particles to overcome a huge potential barrier given by V £ (y ) = \V(y), 
see Figure 15.11 

For simplicity we assume that the total mass f Q u dx dy as well as the volume |D| 
of the physical domain equal 1. Hence, we can again consider the model as a Markov 
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Figure 5.1: The potential V along the reaction path T = [0, 7], 


process with continuous paths t (->• (X t , Y t ) G x T = Q, whose distribution laws can be 
described by densities u(t) G Prob(Q). The Kolmogorov forward equation reads 

it — mnA x u + T e dy(d y u + udyVe), (V x u, d y u+u d y V £ ) ■ v — 0 on dQ. (5.1) 

Clearly, the unique steady state w £ is independent of x and takes the form 

d y. 

Equation (15.111 is a Fokker-Planck equation and, hence, has the Wasserstein gradient 
structure introduced in [JKQ98 ] with 

d £ {u) = jJ \ B (^-)w e dydx and %* e (u,£) = JJ V x £| 2 + y(d y £) 2 )udydx. (5.2) 

For studying the limit e —» 0 we now assume that V G C 2 (T) has exactly two non¬ 
degenerate minimizers as pure states, where V = 0 w.l.o.g, and one global maximum as 
barrier, namely 


w e (x,y ) = w E (y) := — exp (-^K(r/)) with := J exp 


V{2) = 1/(6) = 0, V(y) > 0 on T \ {2, 6}, V"{2) > 0, V"(6) > 0; (5.3) 

1/(5) > V(y) on T \ {5}, V"(5) < 0, (5.4) 


see Figure EHJ (Again, any two points in T could be taken as the pure states, and any 
point in between as barrier.) As a consequence w £ concentrates in the points y = 2 and 
y = 6 in the limit, viz. 


w £ 


w 0 = a 0 5 2 + an<f 6 G Prob(T), a 0 


y/v"(6) , _ yT"( 2 ) 

v ' v ’ 1 " 11 ''' 1 y / V /,, (2)+y / V"(6)' 


(5.5) 


Here the convergence means f T <f>(y)w £ (y) dy —> ao4>(2) + cti0(6) for all cf) G C°(T). 

The important point is now to choose the diffusion constant r £ sufficiently large such 
that the transitions between y = 2 and y = 6 can occur on times of order 1. According 
to Kramer’s rule (see e.g. [ AM*12] 1. this is achieved by choosing rriy > 0 and setting 


t £ := my 


—rwdy, with gives — exp (— ViTd/e) 
w £ (y) e v 2 


From the concentration of w £ in the points {2, 6} we obtain that w £ G Prob(Q) con¬ 
centrates in the sets fl x {2} and fl x {6}, namely 


w £ -A Wq ■= <g) wo in Prob(Q). 
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Recalling that £ e is the relative entropy with respect to w £ , the T-convergence £ £ ^* £o 
appears natural. To be more precise concerning densities and measures, we define 

p , _ f E((co, Ci)) for//= Coda;®^ + Cida;®6 6 , 

’ { oo otherwise, 

where E((co,d)) := / (a b (^)o;o + A B (^)ci!i) dx. (5.6) 

Proposition 5.1 We have £ £ £ 0 in the weak* topology of Prob(Q). 

Proof: The lirninf estimate is established in [AGS05 . Lem. 9.4.3]. 

To construct recovery sequences, we may restrict to the case £o(/x) < oo, since oth¬ 
erwise the lirninf estimate provides the result. Hence, we may assume /i = c 0 da;®^ + 
Ci do;®5 6 and, using a nonnegative, continuous cut-off function x '■ V •-> max{l — \y\, 0}, 
we can define the measures 


p £ = u £ (x,y)dx®dy with u £ (x,y) = c 0 (x)w e (y)p 0>e x(y-2) + Ci(x)w £ (y)pi, e x{y- 6), 
where the normalization constants are given by 

Po,s J T w e (y)x(y~ 2) dy = 1 = P he f r w £ (y)x(y- 6) dy, 

which implies Pj t£ —* 1 /oij > 0 for e —> 0. Then, we easily find £ £ (u e ) —> £o(h)- ■ 

Thus, the limit evolution will be described by the densities Co and ci on for the 
particles being in the pure states y = 2 and y — 6, respectively. In particular, in the limit 
e —y 0 the time that the particles spend along the reaction path away from these points, 
i.e. in T \ {2, 6}, is 0. 

One difficulty in deriving the lirninf estimate for De Giorgi’s dissipation functional 

S> £ (u) \= I %(u(t),u(t)) + %* E (u(t), -DE e (u(t))) dt 
Jo 

is that tR £ is only implicitly defined via the Legendre transform of Jl*. Moreover, we are 
not able to employ the classical Wasserstein gradient flow theory in [AGS05J using the 
Benamou-Brenier formulation, because of the different roles of the diffusion in x with 
mobility mo and the diffusion in y with mobility t £ — » oo. The first step to establish 
the following result follows the idea in [ MaM15a j. where one obtains a lower estimate 
by replacing ‘R £ {u,ii) by the smaller term ( £ e ,u e ) — Jl*(u £ ,^ £ ) and by choosing a suitable 
recovery sequence —> £ 0 f° r the limit passage £ —> 0. Finally, one takes the supremum 
over all £o to recover 0l o as dual of OIq. The second step involves a suitable transformation 
of the reaction variable z = Z e (y) (first introduced in [ AM* 12] ) which allows us to control 
the relative densities v £ := u £ /w £ and the dual potentials £ £ along the reaction path T. 

In the following result we will again describe the limit GS (Prob(Q), £ 0 , 3? 0 ) by a 
reduced GS (Prob(G x {0,1}), E.R), since in the limit every /j G Prob(Q) with finite 
relative entropy satisfies y — Co dx®52(y) + c\ dx®5e(?/) with (co da:, c\ da;) G Prob(h2 x 
{0,1}), see (15TD. 

Theorem 5.2 (From diffusion to reaction-diffusion) The family of gradient systems 
(Prob(Q), £ £ , 0l £ ) defined via (15.21) converges in the EDP sense to the gradient system 














(Prob(Q), £ 0 , 3 ? 0 )j where £ 0 is given in (15. 6 p via E and accordingly 0Z o is given via R, 
which is defined in terms of the dual dissipation potential 

R *( c > 1 i) '■= / n ^(co|V^| 2 + Ci|V^| 2 ) +m r ^]^ c tf*(ri 1 -r) 0 )dx. 

The above result means that the limiting GS is a generalized gradient system defined 
for c = (co,Ci) G ProbfRx {0,1}), where the limiting system is the coupled system of 
linear PDEs given in the form 

c 0 = m n Ac 0 - my ( c 0 /a 0 - Ci/aq), c x = m^Acx + m r (c 0 /a 0 - Cx/af) , 

with Neumann boundary conditions Vc.,- • v = 0. We emphasize that the original GS 
(Prob(Q), £ e , 0l £ ) is the classical GS for the Fokker-Planck equation, while the EDP limit 
provides the generalized gradient structure discussed in Section B.4.31 We observe that for 
e > 0 as well as for e = 0we have the GS that is induced by the large-deviation principle 
discussed in Section 12.41 Thus, we have found another instance of the interchangeability 
of the large-deviation principle and the EDP-limit, as displayed in Figure [ITT] 

p 

Sketch of proof of Theorem ] 5. 21 Since the T-convergence £ e —*■* £o was already established 
in Proposition 15.11 it remains to show the liminf estimate for the dissipation functional 
ZZ £ . More precisely, assume u £ (t) —^ /i{t) = c 0 (t ) dxG)<5 2 (2/) + Ci(i) d^c<g)<5 6 (?/) in Prob(Q) 
for all t G [0, T ] such that sup te r 0 T i £ £ (w £ ) < oo; then, we have to show 

liminf ®Ju e ) > D(c) := / R(c, c) + R*(c,-DE(c))d t. (5.7) 

Jo 


Step 1. Dualization of 0l £ : The first major idea follows [MaM15a] and exploits the 
definition of A e as Legendre transform of 01*. Introducing the functional 

&e{u,£)\= [ (£,w) - ft*(tt,0 + 0i*(u, -D£ e (iz))dt, 

Jo 

we easily see that J$ e {u) can be reconstructed via sup^ £$ £ (u, £). Using the definitions of 
£ £ and 01* we have the explicit form 





0 JQ 


fu - ^|V x f | 2 u - ^-{dyifu 


+ 


mn \V x u\ 




u 


u 


d y dx d t. 


Step 2. Rescaling the reaction-path variable. The second major idea follows [ AM* 12 . 
Sec. 2.1], where no x-direction was present. We define the diffeomorphism Z £ : T —> Z : = 
[0,1] and its inverse Y e = Z~ l : Z —> T via 


2 = Z e (y) := 


my 


Y Jy =0 


^d y and Y e '(z) = — w e (Y e (z)). 

my 


The transformed equilibrium density w £ on Z is 

w £ (z) := w £ (Y £ (z))Yf(z) and satisfies w £ wq := ctcAo + «i^i- (5.8) 
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Indeed, for the latter statement we first use that for all g G C°(Z>) we have the identity 
fzd(z) w e (z) dz = f r g(Z £ (y))w £ (y) d y. Recalling that V has a unique global maximum 
at y = 5, the function Z £ converges uniformly on compact subsets of T \ {5} to the step 
function Z 0 (y) = 0 for y < 5 and Z 0 (y) = 1 for y > 5. With this and (15.5|) we conclude 
f r g(Z e (y))w £ (y) dy —> a 0 g(Z 0 (2 )) + a 1 g(Z 0 (6 )) which is the desired result (15.8p . 

To estimate in the limit e —> 0 we use now the independent variable z = Z e {y) and 
the dependent variables 

v(t,x,z) = — - 7 and C (t,x,z) = £(t,x,Y e (z)). 

w e\le\ z )) 

Introducing the domain Q = x Z we find & e (u, £) = & e (v, () with 


@e(v,C) 



(vw £ — —\V x (\ 2 vw e 


mr 


{d z Q 2 vd 


| V :c n| 2 
2 v 


W e d 


m r (d z v) 2 

~2 V~ 


dz dx df. 


The transformation of 

udy — vw £ dz, d y f = 


to £$s follows easily by using the relations 
d z ( , r e mr 


U'W’ 


and 


(Y’{z)) 2 w e (zY 


Step 3. The T-limit for S$ £ (-,Q: The importance of the new form SS £ is that the 
dependence on £ only occurs in the weighting measure w e . Since w £ concentrates in the 
points z — 0 and 1, the three terms that are multiplied by the weight w £ will converge 
to simple integrals over [0, T] x for the densities Co and c\ respectively. In contrast 
there are two terms not involving w e , but these terms only involve derivatives in the z- 
direction. In particular, they control the smoothness of C and v in z-direction, namely 
v / Tgl 2 ([0,t] x fl; H 1 (2,)) such that vw e indeed has a well-defined limit. With this and 
w £ —^ Wq = a 0 8i + Oi\8\ (cf. (15 .8p ). it is possible to show that for fixed and sufficiently 
smooth Q we have & e {-, () &o(-, () with 


&o(v,0 





(d z 0 2 v) dz 


3=0 


E l /■ • m n IV7 >. |2 , '"-a 

a j[Cj V j cy V j\^xCj\ + ry 


% |V x u 


12 1 


dx di, 


where Vj(t,x) = v(t,x,j) and ( j(t,x) = £( t,x,j ) for j — 0 and j = 1. 

Step 4 . Minimization over the reduction path profile. Note that in the definition of 
SSq, the values of the functions v and (, f° r z G ]0,1[ only occur in the first integrand (with 
factor mr). Hence, one can eliminate the integral by taking the supremum in ( and the 
infimum in v for given boundary values at z = 0 and z — 1. The relevant functional reads 


and Proposition IA.3I provides the following explicit inf-sup formula 


inf sup { N(u, C) | C(0) = Co, C(l) = Ci } u(°) = v 0 , n(l) = v lt v > 0 


= y/v 0 vi c &*(log w - log v 0 ) - y/voVi if*(Ci-Co) =: N((i-(o,v 0 ,vi). 
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Thus, we can reduce to a functional B on v = (v 0 ,vi) and £ = (Co,Ci)> namely 


B(v,C): = 



o Jn L 


m r N( Ci—Co, ^ 0 , ni) 


( . m n 2 m n \V x Vj \ 2 

/ y W ( Q] v j „ Til^COI + 9 


3= o 


dx dt. 


The inf-sup dehuition of N provides the following relation between and B: 

V v with u|[ 0 ,T]xnx{o,i} = v3( with CI[o,T]xnx{o,i} = C : 4)(+ () > B(v, £). (5.9) 

Step 5. Identification of the limits: It now remains to relate the limit functions v to 
the weak limit of the sequence u e . For this, we consider a sequence u £ as in (I3.4cp . i.e. 
u £ — 1 u and E £ (u £ (t)) < C < oo. By the definition of v we have 



o JQ 


u £ (t,x,y)(l)(t,x,y)dydxdt = / v e (t,x,Z e (y))w e (y)<j>(t,x,y)dydxdt. (5.10) 

Jo JQ 


Without loss of generality we assume oo > C > S> £ {u £ ) = sup ^U$ £ {u £ ,tf) = sup<-J^(u £ , £), 
which gives the bound 

|| \/’+1|l 2 ([0,T] x+H 1 (2.)) < C- (5-11) 

This implies Hblder continuity of v(t,x, •) : Z —> R. Moreover, Z £ converges uniformly to 
0 and 1 near y = 2 and y — 6. Hence, we can pass to the limit in (I5.10P and obtain 



cT 

i o J n 

= lim 

£—>-0 


cfit, x)cj)(t, x, 2)+c 2 (f, x)cj)(t, x, 6) dxdt = lim / / u £ (t, x, y)4>(t, x, y) dydxdt 

£ ^° Jo Jq 

v £ (t, x , Z e (y)w e {y)<j>(t, x, y) dydxdt = I / E Vj(t, x)aj(p(t 1 x, 2+4 j) dxdt, 



o jq 


o Jn 


which means c = (cq, +) = (aoUb oiiVi). Using the explicit form of R* and E implies 


[ T (C, C>— R*(c, C)+R*(c, —DE(c)) dt = B(( Cl /«o, c 2 / ai ), C). 
Jo 


(5.12) 


Step 6. The liminf estimate: With these preparations we can now complete the liminf 
estimate. By the construction of v £ and f £ (t,x,y) = C(f,£, Z £ (y)) we obtain the relations 


@e{u e ) > & £ {u £ ,f £ ) = 4(wC), 

where ( is now fixed. For the sequence u £ as given in Step 5, we can further assume that 
v £ —- v in L 2 ([0, T\ xfi; C°(%)), using (15. lip . According to Step 3 the liminf for e —» 0 
yields 

liminf %(u e ) > lim inf 4(w C) > 4(+C) > B((c 0 /a 0 , Ci/cq), C), 

where for the last estimate we have to choose ( according to (15.9p to fit the limit v 
and C — CI[o,t]xUx{o,i}- Nevertheless, the functions (, = (Co? Ci) are still free. Using the 
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characterization (I5.12p and taking the supremum over all £ gives the desired lower bound: 

lim inf @ e (u e ) > sup3§o(v,Q > sup / ((£, c) —R*(c, 0+R*(c, — DE(c))) dt 
£ ^° C C Jo 

(R(c, c) + R*(c, —DE(c))) dt = D(c). 

Thus the desired estimate (15.711 is established, which finishes the proof of Theorem 15.21 ■ 



A Evaluation of some functionals 

Here we give explicit calculations for the functional S occurring in the membrane limit 
and the functional N occurring in the limit of diffusion to reaction. It is surprising that 
both functional are closely related, see (IA.5p . 


A.l Derivation of the potential G(a,uo,ui ) 


We first give the result of the standard case of constant coefficients A and W, which was 
already derived in |MPR141 Prop. 4.4] under the restriction u 0 + M i — 1- For the functional 

fl a 2 + u'(x ) 2 


S (a,u) := 


2 u{x) 


■ da;, 


(A.l) 


we define the value function 

G(a, uq, Mi) := min j 9(«, m) u e H 1 (0,1), m(0) = uq, m(1) = u\, u > 0 
and give a full proof of the derivation of the explicit formula. 

Proposition A.l For all a G M and Uq^Ui we have 

G(a,u 0 ,ui ) = + v / MoMi^'*(logM 1 -logM 0 ), (A. 2) 

where the last term simplifies to G(0,m 0 , m i) — 2 (y / %—yfixf) 2 . Moreover, the unique 
minimizer is given by 

u(x) = {1 —x)uq + xu\ + b(x 2 —x) with b = uo + u± — \J a 2 +4M 0 Mi. (A-3) 

Proof: Since the integrand is convex, there is a unique minimizer u. Denoting the 

integrand by f{u,u') the Euler-Lagrange equations —(d u ’f(u,u')) + f(u,u') = 0 are 


uu 


u 


+ cm 2 = 0. By Noether’s theorem we also have the first integral u!d u ’ f{u,u') 


f(u,u') = ((w') 2 —a 2 )/(2m) = y/2 =const. From (m ') 2 = a 2 + 7 M it is now easy to see 
that all solutions of the Euler-Lagrange equations are parabolas. Using the boundary 
conditions we find m in (IA.3I) . where 7 = 46. 

To evaluate the integral we restrict to the case u'(x ) > 0 on [0,1], which means 
2 mo < \/« 2 +4moMi < 2mi. I 11 the other cases, one can do the calculation on all monotone 
parts in a similar fashion and add the result. We use (1A.3|) and (m ') 2 = a 2 + ju to obtain 


a 2 + u'ix) 2 7 2 f 1 dx 

i 2 u(x) 2 J 0 u(x) 


fUl 


= 2b + a 1 


du 


f UQ 


U\J q ? 2 + 7 m 


= 2b — 2a artanh 


T^Mi + 2 a artanh ^ 1 +^mq. 
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To proceed we first observe a 2 + 7 Uj = (\/a 2 +4 u 0 u\ — 2Mj) 2 , which gives \J a 2 +7 m 0 = 
\/a 2 +4'UoWi — 2m 0 and ^ a 2 + 7 Mi = 2mi — vA 2 +4moMi. Now employing the addition rule 
artanh(x) + artanh(y) = arsinli [fx-\-y)/\J (1— x 2 )(l— y 2 )) and 7 = 46 gives flA.21) . ■ 

In Section |4] we need a more general version with non-constant functions A and W: 

S(-.f') := £Sj + ^((lo g (U/W))') 2 dy. (A.4) 

We will show that the influence of the coefficient functions A and W can be calculated 
from Proposition I A. II bv a suitable rescaling of the layer variable in the form x = X(y). 


Proposition A. 2 We have the following formula: 

G(a, uo , Mi) := min{ 9(a, U) \ U > 0, U(0) = Mo, 17(1) — u\} 

= + A*,/^W log(«)V 

* w 0 wi \A, y u 0 ui J * y wowi y ^yuiwo/J 1 

where w 0 — W(0),wi — W(l), and A* = Harm(AW) = ( f* l/(A(y)W(y))dy) 1 


Proof: We dehne the new independent variable z and a new function v(z) via 


z = z (y ) : = -4* [ 

Jo 


dy 


0 A{r])W{r]) 


and v(Z(y)) = 


u(y) 
w(yy 


where A* = ( 1/(A(rj)W (?])) dy) \ By definition we have Z( 0) = 0 and Z( 1) = 1, and 
the inverse Y of Z maps [ 0 , 1 ] into itself again. Hence, using Z'(y) = A*/(A(y)fP(y)) the 
functional 9 from (IA.4D is transformed into S via 9(cy U) = A* 9(a/A*,u). The result of 
Proposition IA.2I now follows from Proposition IA.1I via 

G(a, Mo, Mi) = min{ 9(a, U) | C7( 0 ) = Mo, 17(1) — u\ } 

= min{ A*9(a/A*, v) | m(0) = U(0)/W(0), u(l) = U(l)/W(l) }. 


Thus, the asserted formula is established. 


A.2 Derivation of the potential N 

We consider the functional 

for functions v > 0. Hence, N is convex in v and concave in £. We are interested in the 
inf-sup for given boundary values, namely 


N(5, Vq, Mi) := inf { sup { N(u, C) | C(1)~C(0) = 5 } u( 0 ) = v 0 , v(l) = v 1} v > 0 


The following result provides an explicit formula in terms of the dual dissipation 
potential c £*. It is based on the following surprising relation between N and 9 from dA.lj) : 


M(h, m) := max { N(m, C) | C(l) — C(0) = } = nun { 9(ct, v) — aS | a G M }. (A.5) 
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The equality = can be checked by elementary calculations, since in both cases we find 

f 1 dz n -i 


v r (z\^ 6 ^ 

M(5, v) = / — - - d z -Harm(v), where Harm(t;) = 


v(z) 


2 v(z) " 2 

Using the strong link (IA.5|) between K and S we show that N can be calculated from G. 

Proposition A.3 We have the relation 

N(5,V o,Ui) = y/V 0 Vi (log Vy- log V 0 ) - y/VQVx <£*(5). 

Proof: Using (I A. 5ft we want to show that N is related to the Legendre transform 

G*(S,v o,ui) := sup agR <5o — G(a,vo,vi) of G from flA.ljl . For this we keep 8 G R. fixed. 

The functional (a,v) i—>• 5(a,v) — is jointly convex, such that it can be minimized 
in any desired order of a and v. Letting V := { v \ v > 0, u(0) = vq, u(1) = v\ } we have 

N(8,Vq,Vi) = inf M(h, v) = inf ( inf S(cr, u)-”AaA 
v&V v£V VaGR 7 J 

= inf ( inf S(a,n) —= G(8,v 0 ,v i) — Sa = —G*(8,u 0 ,ui). 

qEM V v£V / aGM 

Thus, evaluating G* with G from (IA.2[) explicitly gives the desired result. ■ 
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